Non-integrability of density perturbations in the FRW universe 
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We investigate the evolution equation of linear density perturbations in the Friedmann- Robertson- 
Walker universe with matter, radiation and the cosmological constant. The concept of solvability 
f~>) . by quadratures is defined and used to prove that there are no "closed form" solutions except for the 

known Chernin, Heath, Meszaros and simple degenerate ones. The analysis is performed applying 
5^ ■ Kovacic's algorithm. The possibility of the existence of other, more general solutions involving 

a 

special functions is also investigated. 
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I. INTRODUCTION 



This paper is an attempt to apply the methods of differential algebra, and the differential Galois theory in particular, 

OO 

to a problem of cosmology. Although the considerations are mostly mathematical, the problem itself, and its solutions, 

10 ! 

are of rather practical interest. Linear perturbations of the Einstein equations in many forms are investigated due 
O to their direct relation to such practical questions as the formation of galaxies or the CMB inhomogeneities. We, 
| ' however, enter the physical domain only as the source of a theoretical problem, on which we concentrate. 

The result itself is of negative nature, or, in other words, it makes any further investigation of this kind unnecessary. 



All the known solutions are given, together with the conditions for their validity. As no new ones can exist, this closes 
and completes the analysis of the given equation. 

Of course, that is not to say that apart from those special cases nothing can be said about the behaviour of the 
solution. It is only to some extent that differential algebra can make exact the intuitive concepts of "being possible to 
solve" or "expressible in a closed form" . The definition of non-integrability employed here is but one of many which 
were born as classical mechanics evolved. The Liouville's theorem implies that enough first integrals might yield a 
complete solution of a dynamical system, and accordingly many criteria regarding the existence of certain classes of 
first integrals were developed. The first achievements were those of Kovalevskaya and Lyapunov, greatly improved 
only recently by Ziglin, Morales and Ramis. The Galois theory used here, can also be applied to prove non-existence 
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of meromorphic first integrals in more complex systems. Paradoxically here, we are presented with an equation simple 
enough not to allow for the application of those advanced methods. Thus, only a small part of the theory is put into 
practise and explained here. For a complete bibliography and exposition see for example [7]. 

The paper is organised as follows. In section 2 we derive the equation in question in a non-standard but intuitively 
clear way. The next two sections describe the concept of Liouvillian solutions, integrability and give the basic criteria, 
which we proceed to use in section 5. We also investigate the possibility of solving the problem by a combination of 
some special and Liouvillian functions, and give the theorem which is the main result in section 6. Finally, conclusions 
and finishing remarks are given in section 7. 



II. DENSITY PERTURBATION EQUATION 



We will be considering the Friedmann- Robertson- Walker universe given by the metric 

dr 2 



ds 2 



c 2 dt 2 - a{tf 



1 - Kr< 



+ r 2 dtf 



(1) 



where K is the curvature index, dfi 2 the distance element on a two-sphere. The universe will be filled with radiation 
and baryonic matter characterised by their pressures and densities p and p. A non-zero cosmological constant's effect 
will also be considered. 

The Einstein equations for this model give 
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where G, and A are the gravitational and cosmological constants respectively. The conservation equation reduces to 



c 2 p= -3H(p + pc 2 ), 



(3) 



with the dot representing the time derivative. This can be expressed as the following transformation laws for matter 
and radiation respectively: 

PmO v a 



PmO v a / 
Pr = f^£\ 4 
PrO \ a ) 



(4) 



r — 3fr 



PmO and p r0 are the values of the densities for the moment when a = a , which can be chosen as the present day. 
The fluctuation is introduced by the means of the scale factor 



a = a(l + y), 



(5) 



and a is the solution of the original equations (2). As it is a matter perturbation only, we have 

Pr = Pr, 
Pra = Prn (|) 

= p m (l + y)- 3 
= p m (l - 3y). 

Where we use the scaling law (4) and linearise the problem. Substituting this into the second of the equations (2), 
we obtain 



5(1 + y) + 2&y = -^f-~a(l + y)p r - ^f-~a{l + y)p m (l - 3y) + ±Aa(l + y), 
which, after substituting the equation satisfied by the unperturbed a, simplifies to 

y + 2Hy - A*Gp m y = 0, 

with the Hubble "constant" 

H:=i 



(6) 



(7) 



(8) 



This kind of perturbation is of the scalar type, i.e. constructed from a single function y. Here, we take it to 
depend on time only, although in general it could also involve spatial variables. This case could then be thought 
of as the zeroth order coefficient in the expansion of y(t, r, 9, <p) in terms of eigenfunctions of the spatial Laplace 
operator. Another generalisation would be to consider the vector and tensor type perturbations (see [11] for details 
of the decomposition). As it turns out, vector perturbations also admit quite general exact solutions [10]. 

In order to simplify the equation (7) fully, that is, bring it to the linear form with rational coefficients, we choose 
new variables 



— , u := H t, 
ao 



and constant density parameters: 

flxo '■= 
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which allow us to rewrite the main equation as 
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(9) 

(10) 
(11) 

(12) 



x \du J 

where the perturbation y is now considered as a function of x. Finally, using the first of the equations (2), which in 
the new variables reads 



da 

— ) = sIaqx 4 + n K0 x 2 + n m0 x + n rQ =: W(x), 



(13) 



we can eliminate the derivatives with respect to u, and denoting the differentiation with respect to x with a prime, 
we get 

x{Q.K x i + il K0 x 2 + il m0 x + VL r0 )y" + (3f2 A o£ 4 + 2fl K0 x 2 + §0 TOO a; + fl r o)y' - § M m0 y = 0. (14) 

As follows from the definitions, Qko and f^o are of arbitrary signs, while fi r o and fi m o are non-negative. We take 
TO o to be strictly positive, though, because of the nature of the examined perturbations. 

We also introduce the functions p(x) and q(x) related to the above equation in the following form: 

y"(x) + p(x)y'(x) + q(x)y(x) = 0. (15) 

They are 

, x _ 6^ A0 x 4 + 4tt K0 x 2 + 3£l m0 x + 2fl r0 ^ d -3^ m0 

2x(£l A0 2; 4 + flKOX 2 + QmOX + ^ro) ' 2x(tt A0 X i + QkoX 2 + ttmOX + ^ro) ' 

Having obtained the solutions as functions of x it is straightforward to express them as functions of the cosmological 
or conformal time, since x(u) satisfies equation (13), and is therefore expressible in terms of the elliptic functions. 
The exact formulae can be found for example in [8] . 



III. LIOUVILLIAN SOLUTIONS 



Equation (14) is, in general, a Fuchsian one, and the solutions can be found by means of series. In particular, in 
vicinities of the singular points 

y±(x) = (x- x ) a± (l + £ c n{x - x ) n ^j , (17) 

where a± are the characteristic exponents at the point xq. Such solutions are local, and the area of convergence 
around a given point is restricted by the remaining singular points. However, there are some special cases in which 
the solutions can be expressed by means of known special functions, and become global. It is then much easier to 
investigate and understand their behaviour. Thus, we are lead to the natural question of existence of such "simple" 
solutions. In this section we give a short description of a class of functions, which could here be called closed- form 
solutions. 

It is natural to start seeking for the solutions in a class of functions to which the given equation's coefficients belong, 
but such a set proves to be insufficient in most cases. As we are dealing with an equation whose coefficients are rational 
functions, the first choice is the field C(x) - rational functions with complex numbers as the field of constants. Or, in 
the language of differential algebra, (C(x),d) - the above field equipped with a suitable derivation operation, which, 
in this case, coincides with the usual, complex one, and will be denoted by a prime throughout this section. 

This field is much too small, of course, and we are soon forced to extend it, introducing new elements. Just, as 
C(x) is an extension of C obtained by adjoining an indeterminate variable x, our new fields will be <C(x) to which 
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new functions are added. Naturally, the derivation on the extended field must coincide with the subfield's derivation 
when restricted to it (that is essentially what a differential extension is). 

To keep the new functions relatively simple, so that the usual notion of "solvability by quadratures" could still be 
applied, the new elements are restricted to three classes. 

Definition 1 For a differential field extension F C G, an element a G F is: 
primitive over G if a' e G, 
exponential over G if a' /a e G, or 

algebraic over G if P(a) = for some P(x) G G[x] - the ring of polynomials with coefficients in G. 

Definition 2 A field extension is called Liouvillian if it is a result of a finite number of extensions, each the adjunction 
of an algebraic, exponential, or primitive element. 

Some examples of "new" functions appearing in this process are radicals (for algebraic elements), logarithms and 
inverse circular function (for primitives), and trigonometric functions (for exponentials). In short, the new elements 
are expressible as combinations of exponentials, integrals, and algebraic and elementary functions. 
Finally, we are able to formulate what we mean by intcgrability: 

Definition 3 A linear differential equation with coefficients in C(x) is said to posses Liouvillian solutions, or be 
solvable by quadratures, if its solutions belong to a Liouvillian extension ofC(x). 

IV. MONODROMY AND GALOIS GROUPS 

In order to characterise a differential equation we can introduce two groups whose invariant properties are closely 
connected with those of the first integrals of the given equation. 

The first, called the monodromy group M. is associated with the continuation of the local solutions of a linear 
differential equation, along closed paths in the domain where the equation is defined. The group itself is an image the 
fundamental group of that domain, and a subgroup of GL(n, C) for n the order of the equation. The only other fact, 
concerning M, that we will need here is that M. C Q the differential Galois group, which we proceed to describe. 

Q is directly connected with the extension of the base field C(x) to a bigger field F, which contains the solutions 
of the considered equation. Q is defined as the group of automorphisms of F that leave C(x) element-wise fixed. 
Carrying this information regarding the equation, the Galois group enables to test for particular forms of solutions. 
One of the fundamental properties is as follows: 

Lemma 1 Let Q be the differential Galois group of a linear differential equation w"{x) — r{x)w{x), r(x) e C(x). 
Then it is an algebraic subgroup of SL(2,C), and one of the following cases can occur. 
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1. Q is triangulisable. There is a solution of the form exp J ui, where u G C(x). 

2. G is conjugate to a subgroup of 



ci 
c{ x 



ci e C* } u 



c 2 

-C^ 1 



c 2 e C* 



w/iere C* = C \ {0}. There is a solution of the form exp J w, with lo algebraic over C(x) of degree 2. 
3. G is finite. All the solutions are algebraic. 
4- G = SL(2,C). The equation has no Liouvillian solutions. 

V. NON-INTEGRABILITY OF THE EQUATIONS 

We now proceed to the direct analysis of the integrability of equation (14), based mainly on the Kovacic's algorithm 
[1]. First we consider the equation with ^ao = 0, for which the calculation is not too cumbersome, and next, we 
outline the reasoning for the general case. 



A. Models without the cosmological constant 



In order to apply the algorithm, equation (15) is cast into a reduced form 



with 



1 



v"{x) = r(x)v(x), 



3(2kx+1) 2 4kx + 7 

+ 



Ax 2 (kx 2 + x + g) 2 x(kx 2 + x + g) 



Where we have introduced new parameters 



KO 



n 



8 := 



r0 



o ' K ' o 

The reduction itself is performed by means of the following change in the dependent variable 



y(x) = v(x) exp 



p(s)ds 



(18) 



(19) 



(20) 



where the constant xo is arbitrary. This transformation does not change the class of solutions we are interested in, as 
it uses the "admissible" operations only. 

In the general case, there are two distinct roots of the polynomial W(x) in equation (14) and they are both non-zero. 
The degenerate cases will be treated separately. 



Looking for the local solutions around x = 0, which is a regular singular point of the equation with both characteristic 
exponents equal to 5, we find 

vi(x) — x 1 / 2 wi(x), 

(21) 

v 2 {x) = ]n(x)vi(x) + x 1/2 w 2 (x), 

and wi and w 2 are holomorphic at x = 0. Using these, it is now possible to obtain an element of the monodromy 
group, by continuation of the fundamental solution matrix along a small closed path around x = 0. 



(22) 



vi v 2 \ 1-1 -2iri \ I vi v 2 

v[ v' 2 ) \ -1 J \ v[ v' 2 

Since M. C Q, the Galois Group cannot be D\ as that group does not contain non-diagonalisable matrices. Also, 
a triangular, non-diagonal matrix cannot generate a finite group, so that Q itself cannot be finite. We have thus 
excluded cases 2 and 3 of lemma 1. 

Case 1 can also be easily excluded applying the aforementioned algorithm. In order for the solution to be of the 
form 



P(x) exp 



J uj(x)da 



(23) 



with a monic P(x) <E C[x], and lo(x) e C(x), the following equation must be satisfied: 

P"{x) + 2u(x)P'(x) + \uj\x) + uj(x) 2 - r(x)]P(x) = 0, (24) 

for P{x) and u)(x) found according to the algorithm, as follows. First, we define a set of auxiliary quantities: 

o£ = \± 5\/TT46, (25) 

where c € V U {00} - the set of all the finite poles of r(x) including infinity, and b is the coefficient of (x — c)~ 2 in 
the Laurent expansion of r(x) around a given point x = c. Next, for all possible combinations of signs s(oo),s(c) we 
compute the possible degrees of P(x) as 

rf = " s i° c) -E^ (c) J (26) 
cer 

and, for the same combinations of signs s(c), the respective possible uo 



^) = E^- ( 27 ) 



cer 

Inserting the functions just obtained into equation (24) we are left with a system of linear equations determining the 
unknown coefficients of P (x). In this particular case, the only non-negative value of d is 0, so that P{x) — 1, and the 
system is reduced to 

3 



2x(kx 2 + x + q) 



0, (28) 



which cannot hold. The first case of lemma 1 is also excluded, which means that Q is SL(2, C), and the equation has 
no Liouvillian solutions in the general case. 

When we admit a double non-zero root of W(x), or in other words, when 



we have the well known Chernin solution [2] 

y(x) = 2 F 1 (^-iVs,iVs,i;- 



(29) 



(30) 



x + 2Q r0/ 

In the current notation it means that Agn — 1, and that the root itself is xo = —2g. As equation (15) now has three 
regular singular points: 0, x and oo, it becomes a Riemann P-equation. The complete set of solutions is denoted by 



y(x) = P < 



-2q oo 
-iV3 x 
iy/3 1 



> . 



(31) 



Such a P-function is not Liouvillian cither, as can be immediately checked using Kimura's Theorem [4]. 

Another degenerate subcase occurs when g = 0, and k ^ 0, so that becomes a root of W{x) and there is another 
non-zero root x = —-. Like before the solution is a P-function 







y(x) = P 



-- oo 

K 

a; 
\ 1 



(32) 



This time however, it is a Liouvillian solution because we can express the above symbol as two base solutions: 



(x) = X 3/ V^mO + QkoX, 



K0 



y2(x)=x 2 F 1 (l,2,l;- 

where 2 Fi is the Gauss hypergeometric function. 

The next subcase is n — and only two simple roots remain: and —g. We accordingly get 



(33) 



y(x) = P < 



— g oo 

-1 x 



I 3 

u 2 2 



(34) 



As before, this is a Liouvillian function, according to Kimura's theorem, and it can be rewritten as the following two 
independent solutions: 



yi (x) = 1 + — — x, 



2n 



rO 



2/2(z) = y/ n r0 + fl m0 x 2F1 ^2, -i, |; 1 



9. 



rnO 



rO 



(35) 
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which are the solutions discovered by Meszaros [9]. 

The last possibility is that g = n = which implies O m0 = 1> and we simply obtain 

Vi{x) = x, 

(36) 

y 2 (x) = .x 3 / 2 . 

Taking into account that condition (29) makes it impossible for the physical cases (fi TO o > 0) to have zero as a 
triple root, we conclude that the only Liouvillian solutions of equation (14) for fi A0 = 0, appear when at least one of 
the parameters Sl r o or Qko is equal to zero. 



B. Models with the cosmological constant 

The reasoning for f^o ^ is the same as in the previous section, but since the leading coefficient now contains a 
polynomial of the fourth degree, the calculations are somewhat more laborious. 

In general, when W{x) has only simple roots, and they are all non-zero (which means fi r o 7^ 0), we obtain, as 
before, a triangular clement of the monodromy group M.. This leaves us with only the first case of the Kovacic's 
algorithm to check, and the physical requirement of f2 m o > causes the equation to be non-integrable. 

In fact, even if the roots become multiple, but still non-zero, the local solutions in the vicinity of x = do not 
change and, again, we only need to consider the first case of the algorithm. We introduce the roots by the following 
formula 

W{x) = fUo(a; - ei)(x - e 2 )(x - e 3 )(x - 64), 

Taking e 2 = we could rewrite the polynomial as W(x) — $l\o{x — e\){x + el + 2e 2 )(x — e 2 ) 2 , because we must 
have ei + e 2 + e 3 + e± = 0. 

This time, the coefficients b, as defined in the preceding section, is not numeric, but depends on the roots. We now 
have 



1,1/ V3 3ei ei \ l±n 

^ = 2 ± 2V 1 + 4 (4 + ^-^+3^)= : — ■ (3?) 

We see that n must be at least a positive integer, and that initially there are countably many possibilities to check. 
Since the relation between the roots and the densities ft is known, we might obtain an important restriction of the 
form 

n r0 3(n 2 -4) 4 

n A0 " n2 + 12 C2 ' [6 *> 

which means that f^Ao must be negative if n > 2. One can readily check that n = 1,2 give no solution. Further 
calculations with these restrictions on ^ao an d n, determining the coefficients of the appropriate polynomial P(x) 
require us to solve a system of n homogeneous, linear equations. Unfortunately we have been unable to do that for 
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general n, but it is easy to check for the first 50 values that the determinant of the system is not zero, with its modulus 
monotonically increasing with n. Thus we conjecture that there are no values of n for which a solution exists. 

When the double root becomes triple, it also becomes a pole of the third order of the function r{x). This means 
that only the second case of Kovacic's algorithm needs to be considered, but it yields no Liouvillian solutions. 

It is impossible for W(x) to have two double roots or a quadruple root, because as we noted the sum of all roots 
must be zero, and that would lead to f2 m o = 0. 

Letting O r o = changes the multiplicity of x = as a singular point, and the monodromy argument no longer 
holds. We note also that no further increase in the multiplicity of that point is possible as fl m o 7^ 0. 

Assuming first that there arc no multiple roots of W(x), the algorithm immediately yields the following fundamental 
solutions 



yi(x) = x 3/2 v / Qaox 3 + Qkox + SlmO, 

,/ x y/2 (39) 

= 7 {n A0 x3 + n K0 x + nJ dx > 

which holds even if one of the roots becomes double. This solution is also known, and was found by Heath [3]. 
As above, a triple root would mean Q m0 = and is thus not physical. 



VI. ON MORE GENERAL SOLUTIONS 



Although the concept of Liouvillian solutions gives us a simple and applicable formulation of solvability by quadra- 
tures, it is easily seen to be insufficient on itself as a mean to discard an equation as insolvable in general. The fact 
that the Bessel or hypergeometric functions (except for some special cases) are not Liouvillian is the best example of 
that. 

Of course, there are no algorithms for finding more complex solutions, but it is possible to extend the considered 
class somewhat. Following the paper of Bronstein [6], we try to find solutions in the form 

yi (x) = m(x)F 1 [^(x)}, 

(40) 

y 2 {x) = m(x)F 2 [£(x)}, 
where -F\(£) and F 2 (£) are fundamental solutions of a given target equation 

v"(0 = «(0iK0. (4i) 

and m(x) and £(x) are Liouvillian over C(x). The quoted paper presents an algorithmic approach for target equations 
which have an irregular singularity at infinity, and £(x) £ C(x). It also offers a check for a certain class of algebraic 

Upon substituting the form of solutions (40) into equation (18), and using the target equation (41), we are left with 
an expression containing only F(£) and F'(£). Making the coefficients equal to zero gives the fundamental equations 
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to be solved 

3£" (x) 2 - 2?"(x)?(x) + 4 M [C(x)]e'(x) 4 - 4r(x)C(x) 2 = 0, (42) 
mix) = . (43) 

The main theorem of [6] makes it possible to find all rational solutions of equation (42), by bounding the degrees of 
the polynomials involved. 

Theorem 1 Let \\ i Q\ be the square-free decomposition of the denominator of r G C(x). If the order of u(x) at 
infinity ^oc(u) < 2, then any solution £ G C(x) of (42) can be written as £ = P/Q where 

= 11^00(^+2 e CM, (44) 

i 

and P € C[x] is such that either deg(P) ^ deg(Q) + 1 or 

deg(P) = deg(g) + IZTuY ^ 
For the algebraic case, they also provide a possible (but not exhaustive) ansatz of the form 

ax) = p(>(2-~ h)) jjq(*- 2 )( 2 -"oo m) ) (46) 

and a bound for deg(F) of cither 

deg(P) < (2 - ^ 00 (u))(deg(Q) + 2) 

or 

dcg(P) = (2 - ^ 00 (u))deg(g) + 2 - Voo {r). 

Here, we choose to investigate the Bessel and Kummer functions, as the solutions of the target equation. The 
other classical classes of o-Fi and \F\ functions, to which theorem 1 is applicable, are equivalent to these two classes, 
rationally or algebraically. 

In the general case of f^Ao ^ 0, and no multiple roots of W(x), the denominator of r(x) is 16x 2 W(x) 2 , and hence 
its square-free decomposition has only one term Q 2 = xW(x). We take first the Bessel equation, 

F "(0=(^ 1 -^jF(0, (47) 
where e is 1 for the Bessel and -1 for modified Bessel functions. We have Voo(u) — 0, and Voo(r) = 2, so 

Q=n^+2=i, (48) 

i 

and deg(P) < 1. Finally, substituting = c\x + c into equation (42) yields no non-constant solutions. 

Proceeding in the same way, for Kummer functions, algebraic forms of £(x), and possible cases of the roots of W(x), 
we found no new solutions of this more general class. 

Thus, we can finally formulate the main results as the following theorem: 
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Theorem 2 If il m0 > 0, O r0 ^ 0, ftxo <= R> arl ^ ^ao € R (unless W{x) admits a double root, which requires the 
assumption of Qao ^ 0), then the matter density perturbation equation (14) does not posses any solutions Liouvillian 
over C(x), or, in other words, is not solvable by quadratures, except for the following three cases: 

1. Heath's solution (Q, r o = 0)- It is given by (39), for Oao 0, and by (33), when J7ao = 0. 

2. Meszaros's solution (CIko — 0, ^ao = 0). Given by (35). 

3. Matter only (T2ao = ^ro = ^ko = 0, Q m o = lj- Given by (36). 

Moreover, there exist no solutions of the class m(x)F[£(x)], where is a classical special function of the type qF\ 
or iFi, m(x) is Liouvillian over C(x), and £(x) € C(x) or is algebraic of the form (46). 

However, when f^Ao = and W(x) admits a double root, the non- Liouvillian solution of Chernin's, given by (30), 
applies. 

VII. CONCLUSIONS 

The problem studied in this paper is a good example of a possible application of the differential Galois theory in 
practise. The known solution of the perturbation equation were originally discovered in different contexts and over 
the span of a few years, whereas here, one theory allows for full analysis. 

Furthermore, the existence of other closed- form solutions has been ruled out. The one exception is the case of a 
double root of W{x), when the cosmological constant is negative. Although A is usually assumed to be positive now, 
it still remains a viable mathematical possibility. However, we expect that that case will not yield any new solutions, 
as the "manual" check for the first, simplest candidates failed. 

Despite the fact that the analysis of perturbations can be performed numerically, without the need of explicit 
solutions, we feel that a more exact approach is always valuable - providing a better understanding of the complexity 
of the problem, revealing hidden, special solutions, or simply providing a more exact formula for numerical work, ft 
also reflects, in some sense, the fact that general relativity, as the next step in the theory of gravitation, introduces 
non-integrability into the basic equations. This can be seen clearly in a more advanced employment of the Galois 
theory applied to cosmological models themselves, for example in [7] and [5]. 

From the physical point of view, this model is a relatively simple one, and it is already known that the solutions 
found are not on themselves strong enough to explain the density fluctuations levels today. The main problem is the 
fact that they would need to be exponential, whereas they grow at most like a fixed power of the scale factor, and there 
is not enough time for the initial perturbations to increase sufficiently. One has also to be aware that linear instability 
might be dumped when considered in the quadratic regime, so the least to be required of the linear equation's solutions 
is appropriate amplification. Nevertheless, these basic effects, being the first approximation, have physical meaning 
and should be taken into account when constructing more complicated models or explaining observational data. 
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We hope to investigate this field further in the future, as it does not limit the results to negative statements on 
integrability only, as we tried to demonstrate in this paper. 
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